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SUMMARY 


Asymptotic  expansions  of  the  distributions  of  some  classification 
statistics  and  associated  probabilities  of  misclassificatlon  are  con- 
sidered for  a two  population  classification  problem  when  the  population 
distributions  follow  a stationary  Gaussian  process.  Special  cases  have 
been  considered  when  the  population  distributions  follow  a first  order 
autoregressive  process  and,  in  particular,  the  probabilities  of  mis- 
clas8iflcation  is  studied  as  a function  of  the  measure  of  dependence 
between  the  two  populations. 
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1.  Introduction:  Let  u>q  be  an  experimental  unit  which  la  a random 

outcome  from  a population  tt.  It  la  known  that  tt  la  Identical  to  one 
of  the  two  apedfled  populations  tt^  and  V where  tt^  and  denote  the 
same  population  it*  at  two  different  pointa  of  time  t^  and  > reapectively. 
Let  XQ  - X(o)q)  be  a pXl  vector  of  meaaurement  on  the  unit  uj^.  The  problem 
la  to  Identify  tt  with  one  of  tt^  and  tt^  based  on  Xq  and  the  knowledge  of 
the  distributions  of  XQ  In  tt^  and  tt^,  which  are  not  completely  known. 
Information  about  these  distributions  Is  obtained  from  a sample  of  N 
units  0)^, . . . ,ui^  (called  training  sample)  from  ti*  with  X^  as  the 
X-observatlon  on  the  unit  t»a  observed  at  time  or-1, . . . ,K;  1-1,2. 

Let  X^  denote  the  X-observatlon  at  time  t.  We  shall  assume  that 

(1.1)  Xt  - mt  + Ufc, 

where  Ut  follows  a stationary  Gaussian  process,  and  In  particular. 


vhere  E is  a nonsingular  matrix  and  T la  a pXp  symmetric  matrix  (see 
Anderson  (1971)).  Then  («■!»•  «.»H)  are  i.i.d.  and 


E T 


(1.3)  , C ) ~ *2,  !”(  )•(  )] 

ha  * * J 


where 


(1.4) 


*1  " 


1-1,2. 


H I - 
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A special  case  of  (1.2)  will  also  be  Created  when  Ufc  follows  a first  order 
autoregressive  process*  i.e.* 


(1.5) 


°t  ■ XBt-i  + v 


t-0,+1,+2,-"  , 


where  Ixl  < 1 and  e 's  are  i.i.d.  N (0*A).  Then  we  can  take  (see  Anderson 

c p 


(1971)) 


(1.6) 


T - PE,  |p|  < 1. 


Let  denote  the  hypothesis  that  u)q  is  from  tt^  (i«l,2).  When  p^ 
and  Z are  known  the  form  of  a likelihood-ratio  rule  (see  Anderson  (1958)) 
is  given  by  the  following:  Accept  iff 


(1.7) 


W*  = ((X0-n2;S))  - ((Xq-p^E))  > k. 


where  k is  a constant*  and  for  a pXl  vector  Y and  a pXp  nonsingular 


matrix  B 


(1.8) 


((Y;B))  - Y'B_1Y. 


When  some  of  the  parameters  p^,  p-2  and  E in  (1.7)  are  unknown  we  replace 
them  by  their  respective  estimates*  chosen  suitably*  based  on  the  samples 
(X;aX'a) , a*l,... *N.  Such  rules  will  be  called  plug-in  likelihood-ratio 
(PLR)  rules.  The  notion  of  FLR  rules  was  first  introduced  by  Wald  (1944). 

In  this  paper  we  assume  that  ^ and  are  not  known  and  we  consider 
several  cases  depending  on  the  available  knowledge  about  E and  T. 


The  estimates  of  p^  and  p>2  will  be  taken,  respectively,  as  X^  and  X^, 


where 


» 
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N 

(1.9)  X.  - Z X.  , 1-1,2. 

a-1 


Define 

N 

(1.10)  Stj  - S^(Xla-Xt) (Xja”Xj) ' ; 1, j-1,2. 

The  plug- In  version  of  W*  Is  given  In  general,  by 

(1.11)  ^ - (X1-X2)'B"1(X0  - -i^+xp), 
where  the  matrix  B Is  defined  below. 

Case  (a):  B - Z,  when  E Is  known. 

Case  (b):  When  Z and  t are  unknown, 

(1.12)  B - (S11+S22)/(2N-2). 

Case  (c):  When  Z Is  unknown  and  T - pZ  with  known  p, 

(1.13)  B - (S11-pS12-pS21+S22)/(l-p2)(2N-2). 

Case  (d):  When  Z Is  unknown  and  T - pZ  with  unknown  p , 

(1.14)  B - (Su+S22)/(2N-2). 

Note  that  In  each  of  the  above  cases  B Is  an  unbiased  consistent 

estimate  of  Z.  In  naaes  (b)  and  (c),  B Is  an  asymptotically  efficient 

estimate  of  Z;  however,  an  asymptotically  efficient  estimate  of  Z In 
case  (d)  la  difficult  to  obtain  (when  p > 1) . The  limiting  distribution  of 
Wg  a*  N -»  • la  normal  with  variance  a and  mean  If  Xq  is  from  tt^  and  mean 
- If  XQ  Is  from  ^2» 


amaas 


We  shall  assume  a > 0 


In  this  paper  we  derive  asymptotic  expansions  of  the  distributions 

of  % " iV'V72  “d  (Vib  + jvV/2  When  ftX0  " “l  and  eX0  " V 

respectively,  in  powers  of  1/N,  where 

(1.16)  % m ((X^X^B)). 

Anderson  (1973)  obtained  the  asymptotic  expansions  of  these  distri- 
butions when  T ■ 0 and  B is  given  by  (1.11).  We  shall  closely  follow 
Anderson' 8 method;  however  we  have  to  modify  and  extend  Anderson's  results 
because  and  X2  are  correlated  and  the  presence  of  t leads  to  a different 
estimate  of  Z. 

We  shall  also  derive  asymptotic  expansions  of  the  distributions  of 
1 1/2  1 1/2 

(Wjj  - •gx)/a  and  (W^  + ^Jf)/a  when  SXQ  » and  &XQ  *»  y-2»  respectively. 
This  would  extend  the  results  of  Okamoto  (1963)  who  derived  these 
asymptotic  distributions  when  T ■ 0 and  the  usual  estimate  is  taken  for  Z. 
For  deriving  these  results  we  shall  again  follow  Anderson's  method  which  is 
much  simpler  than  Okamoto 's  method. 

Finally,  asymptotic  expansions  for  P (W^  < 0)  will  be  obtained  under 
both  the  hypotheses  and  Hj. 

In  classification  theory,  usually  the  training  samples  from  tt^  and  tt^ 
are  drawn  independently,  i.e.,  on  different  sets  of  units.  This  is  the 
first  paper  where  the  classification  problem  is  treated  with  a dependent 
training  sample  which  occurs  quite  often  in  practice.  The  results  in  this 
paper  also  reveal  the  influence  of  the  covariance  matrix  t on  the  distri- 
butions and  the  probabilities  of  misclasslf ication. 


There  exists  a nonsingular  pXp  matrix  L such  that 


(2.1) 


L£L'  - Ip, 


LTL'  ■ D - diag(P1,...>P  )• 


It  can  be  seen  that  and  0^  (in  all  the  cases  (a)~(d))  are  invariant 


under  the  transformations 


(2.2) 


X0  - L(X0-^1) 
Xi  - L(X1-^1) 

x2  - 

su  - 1 V ! 


Hence,  without  loss  of  generality,  we  shall  assume  that 


(2.3) 


Z - I , t » D - diag(p P), 

D P 


1*1  - 0,  |*2  - -6, 


where 


(2.4) 


6'6  - a. 


As  in  Anderson  (1973),  we  define  Y,  Z and  V as  follows. 


(2.5) 


We,  Y/n1/2,  Xx  - Z/n1/2, 


B - I + V/n1/2, 
P 


n - 2(N-1). 


Y . r r 0 \ r2(n/N)(I  -D)  (n/N)(I  -D)y 

«*••>  c2)~»2p[C0>^„W(i;d)  ^ m 
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f <, 

t 


It 


f. 


t 

r 


i 


Let  Y - (Yj Y )',  Z - (Zlt...,Zp)',  V - [V±J] . Let  Jn  be  the 


subset  of  the  senple  space  defined  by 

,1/2 


(2.7)  Jn  - £ |Y±|  < 4<log  n)l/Z,  \zt\  < 2<log  n)l/2, 


|V±J  | < 2 log  n;  l.J-1,2 p}  . 


The  proof  of  the  following  lease  is  given  in  the  appendix. 


2.1: 


-2, 


P(Jn>  - 1 - o(n  ‘). 


Now 


(2.8) 


Ojj  - ((6«/n1/2;  Ip+V/nl/2))t 


and  let 


(2.9)  Gy,  - (8fY/n1/2)'(Ip+V/n1/2)“1(Z/nl/2), 


(2.10) 


62N  " (Ip+V/n1/2)2)). 


Assume  now  £Xq  « 0 (i.e.,  under  H^).  Then 

(2.11) 


'■<“»  < “>  - e,K«fli,1/2'*<!ui>0ai'1/2l- 


where  i is  the  c.d.f.  of  N(0,1).  Now 


(2.12) 


I**1  (oS|1/2'*€UI)/G2H1/2,*<V 


< e x(JnC)  - o(n"2). 


where  x stands  for  the  indicator  function  of  a set. 
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Following  Anderson  (1973),  it  can  be  shown,  for  sufficiently  large 


n and  Y,  Z,  V in  Jq,  that 


(2.13) 


SlCuO^-M^/G^172]  - «(u)  + n-1/2cp(u)C(Z,V) 


+ n"Vu)[D(Y,Z,V)  - (u/2)C2(Z,V)]  + Rq, 


where  <p(u)  ■ d#(u)/du. 


(2.14) 


(2.15) 


C(Z,V)  - (u/2AZ)6'V6  + (6'Z)/A, 


D(Y ,Z,V)  - A“1(Y'Z-6'VZ) 


+ A^u^'VY-fi'V2*) 


- A"3(6'Y6,Z-6'Z6,V6) 


+ uA"4[(7/8)(6'v6)2-6'Y6'V6], 


A - a . 

Rq  is  the  sum  of  three  terms  Rln,  R2n  and  R3n>  where  R1q  is  a homogeneous 

-3/2 

polynomial  (not  depending  on  n)  of  degree  3 in  Y,  Z and  V multiplied  by  n , 

R2n  is  a homogeneous  polynomial  (not  depending  on  n)  of  degree  4 in  Y,  Z and  V 

-2  5 5/2 

multiplied  by  n , and  R^n  is  O((log  n)  /n  ). 

Next,  we  follow  Anderson's  method,  although  we  clarify  his  results 

and  fill  in  some  gaps.  First  note  that  the  expectation  of  a polynomial  in 

Y^'s  and  Zj's  is  0(1),  and 


(2.16) 


0(n"1/2),  6Vjj2k  ■ 0(1) 
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for  any  positive  Integer  k.  We  shall  now  obtain  the  expectation  of 

the  right-hand  side  of  (2.13)  multiplied  by  x(J  )• 

n 

£[6'Zx(Jn)]  “ 0 since  tZ  - 0 and  Jr  is  symmetric  about  0 in  Z^  and 
the  density  of  Z ^ is  symmetric  about  0.  Moreover,  fi[6'V6x(J  )]  - 0(n"3^2), 


since  tV  “ 0, 


|e(6*v6)x(j  )|  - |«'evx(J  G)fi|, 


(2.17) 


fit |v±J |x( |v±j l > 2 log  n)] 


< [fiVij2]1/2P[X(|viJ|  > 2 log  n)]  - 0(n”3/2) 


by  taking  k ■ 3 in  (A. 3).  Thus 


(2.18) 


fi[C(Z,V)X(Jn)]  - 0(n-3/2). 


Next,  by  Cauchy-Schwarz'  inequality 


|fi[D(Y,Z,V)-(u/2)C2(Z,V)]  - e[D(Y,Z,V)-(u/2)C2(Z,V)]x(J  )| 


(2.19) 


< [fi{D(Y,Z,V)-(u/2)C2(Z,V)}2]1/2[P(JnC)]1/2  -o^'1). 


Similarly 


(2.20) 


leRln-£RlnX<Jn>l  - °<a"2>‘ 


Since  the  third-order  moments  of  the  elements  of  Y,  Z and  V are  either  0 
or  0(n  *^2) , fiR^n  ■ 0(n  *^2) . Note  now 


leR2aX<,,n> I ^lR2nl  “ 0(1)* 


I 
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(2.22)  leR3nX(Jn)l  " 0((log  n>5/n5/2>  “ °<n"2>* 


Combining  the  above  results,  we  get 


(2.23) 


«[(>V/24€U])G2)-1/2l 

- §(u)  + n_1cp(u)e[D(Y,Z,V)-(u/2)C2(Z,V)]  + 0(n'2), 


The  following  theorem  is  thus  proved. 


Theorem  2.1;  When  obtains 


$(u)  + n”1cp(u)fi[D(Y,Z,V)“(u/2)C2(Z,V)]  + 0(n‘2). 


Note  now 


• -h.>  V/2  ^ 


(2.24) 


p2[(^+&,>V/2£-“>- 


where  denotes  the  probability  under  H^,  i=l , 2 . Thus  we  get  the  follow- 
ing corollary. 


i 

I 

f ' 


Corollary  2.1.1:  When  Hj  obtains 

Pt(wR+^,)/atIl/2i>‘l 


- §(u)  - n”1cp(u)e[D(Y,Z,V)-(u/2)C2(Z,V)]  + 0(n-2). 


m 


The  above  two  results  are  similar  to  those  of  Anderson's  (1973)  except 
that  we  needed  a modified  proof  to  treat  our  case.  The  presence  of  T now 


I 

i 


2 

leads  to  new  problems  regarding  the  evaluation  of  £D(Y,Z,V)  and  £C  (Z,V). 
We  shall  consider  this  now  for  different  cases. 

Case  (a):  B - 1^.  In  this  case  V - 0 and  |v^j|  < 2 log  n for  all 

i,j  trivially  holds. 

£D(Y,Z,V)  - A_1£(Y'Z)  - A"3C(6'Y6'Z) 

(2-25)  p 

« A“1(n/N)  I (1-p.)  - A‘3(n/N)6'(I  -0)6. 
l-l  1 p 


In  terms  of  the  original  parameters, 

SD(Y ,Z,V)  - A_1(n/N)tr(I  -Z^t) 

P 


(2.26) 


- A_3(n/N)[A2-(^1-p2)'Z_1TE’1(n1-^2)]. 


For  the  special  case  r - pi, 

(2.27)  £d(Y ,Z,V)  - A_1(n/N)(p-l)(l-p). 

Next 

(2.28)  8C2(Z,V)  - A‘2e(6'Z)2  - (n/N). 

For  what  follows,  (n/N)  Is  replaced  by  2.  Theorem  2.1  now  yields  the 
following: 

Theorem  2.2:  When  B - Z and  obtains 

PI(WN  ‘ ^ u]  " *(u)  + n“19(tt)t2A“1(p-l)  - 2A‘1tr(Z"1T) 
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p . 


iZ  ■* 

«■  - 

i • 


l- 


V 

t 


m 

!■  ' 

*•  I 
{' 


*s 


which  reduces  to 


*(u)  + n”1cp(u)[2A~1(p-1)(1-P)  - u]  + 0(n-2), 


when  t - pi. 


Case  (b):  B - (Sn+S22)/(2N-2) . 

Note  that 


e[D(Y,Z,V)-(u/2)C  (Z,V)] 


2A-1tr(Ip-D)  - 2A"36’(Ip-D)6  - u - uA”  8(6 'V  6) 


(2.29) 


+ [(7/8)uA"4-(u3/8A4)]8(6'v6)2  + 0(n"L). 


In  this  case 

(2.30) 


V - (4n)"1/2[B14B2-2nI  ], 


where 


(2.31)  Bt  - 8U  + S22  + S12  + S21  ~ Wpt2(Ip+D),n/2] , 

(2.32)  B2  - Sn  + S22  - S12  - S21  ~ Wp[2(Ip-D),n/2], 


and  Bx  and  B2  are  Independently  distributed.  Hence 


(2.33)  ev2  - (4n)"1[8(B1-eB1)2  + 8(B2-8B2)2] . 

(2.34)  8(6 *V6)2  - (4n)“1[Var(6'B16)  + Var(6'B26)] 


He  shall  reduce  the  above  expressions  using  the  following  1« 


Lemma  2.2;  Let  A be  a pXp  random  matrix  distributed  as  Wp(A,m),  where  A 
Is  a nonsingular  diagonal  matrix.  Then 

(a)  £(A-&A)2  ■ m[A2  + A(tr  A)], 


and  for  any  vector  6 0 


(b)  Var(6'A6)  = 2m(6'A6)2. 


The  proof  of  this  lenma  Is  omitted. 


ev2  - (4n)-1[4(I  +D)2+4(I  +D)[tr(I  +D)} 
P P P 


+ 4(Ip-D)z  + 4(Ip-D){tr(Ip-D)}](n/2) 


(2.35) 


f(pfl)Ip  + (tr  D)D  + D*], 


(2.36)  e(6'V26)  - (p+l)A2  + 6'D26  + 6'Dfi(tr  D), 


(2.37) 


C(6'V6)2  - 2[A4  + (fi'Dfi)2] . 


From  Theorem  2.1  and  the  above  results  we  get  the  following  theorem. 
Theorem  2.3:  When  B • (Sn+S22)/(2N-2)  and  ^ obtains 

Pl0fc  “ 'K)/V/2  ^ ul  “ #(u)  + a”  V(u) 1 ^tr (I  -£~  1t ) - u 

- 2A"3(n1-»»2)'(I*1-r'1Tr“l)(ji1-n2)  - u(pfl) 

- «A'2o»l-tfc2),E‘lTrlTr'1o*1-^2) 


- 
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- uA'2(n1-n2)'r1Tl-1(n1-n2)tr(TS“1)  + ((7/4)u  - u3/4) 

+ [(7/4)uA"4  - (u3/4  44)][(lji1-^2)'Z-1TE‘1(n1-ji2)]2  + 0(n“2), 
which  reduces  to 

*(u)  + n"Vu)[(2M)(p-l)(l-p)  - u{p  + 1/4  + (p-3/4)p2} 

- (u3/4)(l+p2)]  + 0(n“2), 

when  T ■ pE. 

Case  (c):  B - (s11-PS12-pS21+S22)/(l-p2)(2N-2) , and  D - plp. 

In  this  case  2(N-1)B  ~ W (I  ,n),  and 

P P 

V - n1/2[B-e(B)]. 

Hence,  from  Lemma  2.2 

(2.38)  6(6 'V26)  - (pfl)A2, 

(2.39)  6(6 'V6)2  - 2 A4. 

Theorem  2.1  now  yields  the  following. 

Theorem  2.4;  When  T - pZ,  B - (Sn-pS12-pS21+S22)/(l-p2)(2N-2)  and  ^ obtains 

PK"n-iV/aB1/2i“l 

- *(u)  + n"1tp(u)[24“1(l-p)(p-l)  - u(pfl/4)  - u3/4]  + 0(n“2). 


r:mamm mm 


First  note  chat  under 

(3.1)  Pl(^  - -ia)/al/2  < u]  - e*[(w  + fy3n)Gm~l/2} , 
where 

(3.2)  w « uA  + LZ/2, 

(3.3)  G - (Xj+xp'B"1^-^), 

and  G0  is  given  by  (2.10).  Define  U « (U. U )'  by 

zn  l p 

(3.4)  (X1+X2)  - U/n1/2  - 6, 
and  Y and  V by  (2.5).  Then 


<->  CwO-c* 


0.  .2(n/N)(Ip-D) 


2(n/N)(Ip+D) 


Redefine  J by  , where 
n n 


(3.6) 


J*  - {|Yj  < 4 (log  n)1/2,  |U±|  < 4(log  n)1/2, 
|vij|  < 2 log  n;  i»jal>...»p} 


Then  it  can  be  shown  as  in  Lemma  2.1  that 


(3.7)  P(jJ)  « 1 - o(n"2). 

Replacing  Jq  by  J*  and  proceeding  as  in  Section  2 it  can  be  shown  that 


under  H. 


(3.8) 


PIG^  - ■|of)/a1/2  < u] 

- *<u)  + n'l<p(u)e[D*(Y,U,V)  - (u/2)C*2(Y,U,V)]  + 0(n"2), 

where 

D*(Y,U,V)  - (2A)"1[U,Y+6,VY-a,V2fi-U'V61 
- (2A3)‘1[(6'Y-6,V6)(U,6+6'V6-6'Y)] 

(3.9)  -u[(2A2)_1(36'V26+Y'Y-46'VY)  - 3(2A4)_1(6'Y-6'V8)2] , 

(3.10)  C*(Y,U,V)  - (2A)'l(U'6+a'V6-6,Y)  - uA’2/6'Y-6'v6). 

Now,  in  terns  of  the  original  parameters 

8D*(Y,0,V)  - e(6'V6)2[(2A3)"1  + 3u(2A4)-1] 

- e(6'V26)[(2A)*1+3u(2A2)“1] 

+ (2/A3)(h1-u2)  ' (Z"l-ST1«rl)(|fc1-|fc2) 

- (2u/A2) tr (I  -E” XT ) 

P 

(3.11)  + (6u/A4)  (|*1-|a2) ' (Jfl-STlT2f  + 0(l/n), 

CC*2(Y,U,V)  - 2 + 4(u2/A44«M3)(n1-ji2)'(r1-E’1Tr"1)(n1-n2) 

(3.12)  + [ (1/2A)2  + (u2/A4)  + <u/A3)]ft(6'Vfi)2  + 0(l/n). 

Case  (a):  B - 1 . 

' P 

Hare  V ■ 0 and  we  gat  the  following  theorem  from  (3.8)  to  (3.12). 


i 


Theorem  3.1:  When  E le  known,  B ■ £,  end  obtains 


Pl%  - < uj  - «(u)  + n"Vu)I2A“',((i1-n2)'(£*1-£*1T£"1)(t*1-ii2) 


u - 2uA"2tr(Ip-E"1T) 
2.*3» ». 


2(ujh’-*+u  ‘A"j)0fc1-ti2)  • or^r^r1)  oi^) 


+ 6uA~'t(h1-^2)'(£‘1-£"1t£"1)(^1-h2)]  + 0(n‘z), 


which  reduces  to 


when  T ■ p£. 


*(u)  + n'1tp(u)[2(l-p)&~1-u-2(p-3)(l-p)uA~2 


- 2(l-p)(uV2+uV1)]  + 0(n"2), 


Csse  (b):  £ end  T ere  unknown,  and  B is  given  by  (1.11). 

1 1/2 

In  this  case,  P[(W^  - -jcO/a  < u]  under  Hj  can  be  obtained  from  (3.8), 

| 

(3.11),  (3.12),  (2.36)  and  (2.37);  due  to  Its  lengthy  form  the  explicit  ex- 
pression is  omitted. 

(c):  E is  unknown,  t ■ pE  with  known  p.  and  B is  given  by  (1.12 


Prom  (3.8),  (3.11),  (3.12),  (2.38)  and  (2.39)  we  get  the  following. 


Theorem  3.2:  When  H1  obtains  and  B is  given  by  (1.12) 

- ^«)/of1/2  < u]  - l(u)  + n‘19(«)t2(l-p)6“l-u(3p-l)/2 


- A(p-l)/2  - 2(p-3)(l-p)uA~2  - 2(l-p)(uV2+uVl) 


u(A2/4  + u2+uh)]  + 0(n“2) 
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Similar  results  can  be  obtained  when  H2  obtains  using  the  following 
relation: 

(3.13)  P1[  - lor)/®172  < u]  - P2[  (^  + |®)/®1/2  > -u] , 

where  Is  the  probability  under  (1-1,2). 

4.  Asymptotic  Expansions  of  Probabilities  of  Mlsclasslflcatlon. 

Let  us  consider  the  rule  which  classifies  XQ  Into  iff 

(4.1)  > °, 

where  W^  Is  given  by  (1.11).  Then  the  probability  of  the  mlsclasslflcatlon 
(PMC)  for  this  rule  under  Is  given  by 

(4.2)  PL[  (^  - |a) /a172  < -A/2) , 

and  Its  asymptotic  expansion  can  be  obtained  from  (3.8)  with  u - A/2.  It 
follows  from  (3.13)  that  the  PMC  for  this  rule  under  H2  Is  the  same  as  (4.2). 
Using  results  In  Section  3 we  shall  now  state  the  asymptotic  expansions  of 

(4.2)  for  different  cases  under  consideration. 

When  £ Is  known  and  B Is  taken  as  E,  (4.2)  reduces  to 

♦(-A/2)  + n“1<p(A/2)[A/2  + A-1tr(I  -£-1t) 

(4.3)  P 

“(^r^2) ' (s’1“2"1tZ“1)  (H1->*2){A“3+(4A)“1}  ] + 0(n"2). 

In  particular,  when  T - pE,  (4.3)  reduces  to 

#(-A/2)  + n"VA/2)[A(l+p)/4  + A"l(p-l)(l-p)]  + 0(n"2). 


(4.4) 


When  t ■ pZ  with  known  p but  unknown  I,  end  B is  given  by 
(1.13),  (4.2)  becomes 


(4.5)  *(-4/2)  + n-1<p(A/2)[A(p-l)/4  + A(l+p)/4  + (p-l)(l-p)A“l]  + 0(n"2). 

When  t ■ pZ  with  both  p end  Z unknown,  end  B is  given  by  (1.15), 
(4.2)  reduces  to 


(4.6) 


I (-4/2)  + n-1cp(4/2)[4(l+p)/4  + A(l+p2)(p-l)/4 
+ (p-l)(l-p)4”1J  + 0(n~2). 


Remarks : 

(i)  Assume  p > 1.  Then  the  PMC  of  the  rule  given  by  (4.1)  can  be 

studied  as  a function  of  p,  to  the  order  of  approximations  indicated, 
when  t - pZ.  It  follows  from  (4.4)  and  (4.5)  that  the  PMC’s  in  these  two 
cases  increase  (or  decrease)  with  p when  A2  > (or  <)  4(p-l).  It  follows 
from  (4.6)  that  the  PMC  in  this  case  increases  (or  decreases)  with  p 
according  as  p > (or  <)  [2A-2  - (2p-2)-1l. 

(ii)  When  p ■ 1,  the  PMC's  in  all  the  cases  are  Increasing  functions 
of  p.  However,  it  is  shown  in  Bandyopadhyay  (1974)  that  the  exact  PMC 
decreases  as  p decreases  when  4 is  large  and  it  increases  as  p decreases 
when  4 is  small. 

(iii)  When  T ■ 0,  the  results  in  Section  2 yield  Anderson's  results 
(1973)  and  the  results  in  Sections  3 and  4 yield  Okamoto's  results  (1963). 


APPENDIX 


Proof  of  Lemma  2.1: 

P[|Yt|  > 4 (log  n)l/2]  - (|)1/2  f e'1/2x2dx 

n 

where 

fcn  " 4<1o8  n)1/2[2(n/N)(l-p1)]_1/2. 

By  Mill's  retlo  Inequality 

P[|Yj  > 4(log  n)1/2]  < (|)1/2  e_1/2tn2 

Now 

-t  2/2  L 

e n t_1 
n 

-8  log  n/(n-N)(l-p  ) 

" e [4(log  n)i/Z{2(n/N)(l-Pi)} 

8N/n(l-p  ) 

- [n  (log  n)  ,Z]"l[2(l-Pl)(n/N)]1/2(l/4) 

< [n2(log  n)1/2]‘1(2)_1/2  - o(n“2) . 

Hence 

(A.l)  P[|Yj  > 4(log  n)1/2]  < o(n*2). 

Similarly 

(A-2)  PI |Z±|  > 2 (log  n)1/2]  < o(n“2). 
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From  (1.12)  and  (1.13)  it  follows  that  B * + B2> 

(N-l)B.  ~ W [I  ,N-1],  (k»l,2)  but  B.  and  B,  are  not  necessarily 
ic  p p l Z 

Independent.  Let  Bfc  ■ ((B^^)) * Now 


> 2 log  n 


* (/n/2)[Bul  + B21i  - 2]  > 2 log  n 

* (n/2)[Bllt  + B2ti]  > ZAi  log  n + n. 


Hence,  for  0 > 0 


> 2 lo,  a]  < lo«  n + n/J) 


(e/2)(»ul  + «211)(n-i) 


< #-®(/n  log  n + n/2) 

r OOJ-DB...  0(N-1)B 
• E ^e  111  + e j(l/2) 

. .-»(/»  lo,  o + o/2)(l.M)-n/4>  0 < # < 


(A.  3) 


0 ■ k//n  . 


Fix  k and  let  n be  sufficiently  large  such  that  k/fn  - 0 < 1/2.  Then 


P[V..  > 2 log  n]  < #-Wi/2#-k  log  “(1.2k//n  )"n/4. 


But,  since  n/4  ■ (v/n/2k)k/n /2,  we  have 


PIV^  > 2 log  n]  < Constant  e"  °8  n ■ 0(n”k) , 


Similarly 


P[-V..  > 2 log  n]  < 0<n"K), 


Hence 


(A. 4)  P[|V..|  > 2 log  n]  < 0(n"k), 


Now  let  i + j. 


V^j  > 2 log  n 


» (n/4)(B11j  + B21j)  > A log  n. 


Hence 


PIV^  > 2 log  n]  < e-0^11  log  n E e 


(n/4)0(Buj  + B2ij) 


a/rnlnff  . r (n/2)9B...  (n/2)0B...n 

< e’9^1  108  n E [e  liJ  + e 2ijj(l/2) 

- e-0^  108  n(i-02)n/\  o < e < 1/2. 


Again,  for  fixed  k and  for  9 as  in  (A. 3)  and  for  sufficiently  large  n 
such  that  9 » k//n  < 1/2,  we  have 


P[V£j  > 2 log  n]  < e"k  log  n(l-k2/n)-n/4  - 0(n_k). 


Similarly, 


Pt-V^  > 2 log  n]  < 0(n  ), 


Hence , 
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